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Abstract 


The aerodynamic characteristics of airfoils with several flap con- 
figurations have been studied theoretically and experimentally in an 
environment that simulates a wing immersed in the downwash of a hovering 
rotor. Special techniques have been developed for correcting and validat- 
ing the wind-tunnel data for large blockage effects, and the test results 
have been used to evaluate two modem computational aerodynamics codes. 

The combined computed and measured results show that improved flap and 
leading-edge configurations can be designed which will achieve large reduc- 
tions in the downloads of tilt-rotor aircraft, and thereby improve their 
hover efficiency. 


I. Introduction 


The impingement of the wake of a lifting rotor on a horizontal sur- 
face, such as a wing, fuselage, or control surface, degrades the lifting 
capabilities of the aircraft in hover and lou -speed flight. This "down- 
load" or vertical drag phenomenon, 
is particularly important for tilt- 
rotor type configurations, since 
both the downwash velocities of the 
rotors and the affected wing area 
are larger than for conventional 
helicopters. For example, the 
estimated download penalty in hover 
for the XV-15 aircraft (Fig. 1) 
varies between approxirately 5% and 
15% of the gross weight of the air- 
craft, depending on operating con- 
ditions and the setting of the 
wing flaps. 

One practical, operational 
aspect that illustrates the com- 
plexity of the three-dimensional, 
rotational, separated-flow phe- 
nomena is contained in Fig, 2. This figure shows that the minimum download 
does not occur when the flaps are fully deflected (chat is, when the mini- 
mum wing area is exposed to the rotor downwash) but rather when the flaps 
are deflected approximately 60° . As explained in Section VI, this curious 
behavior now appears to be caused by flow separation on the upper surface 



Fig. 1. The XV-15 Tilt-Rotor Air- 
craft in hovering flight, wing flaps 
fully deflected. 
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Fig. 2. The effect of flap deflec- 
tion on hover performance (Ref. [1]). 
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of the flaps and flaperons. How- 
ever, current engineering predic- 
tion techniques give no clue to 
the mechanism responsible for the 
results shown in Fig. 2. In fact, 
they do not even predict the over- 
all effects of the rotor-wing 
interference adequately, nor do 
they provide reliable design guide- 
lines for reducing the rotor- 
induced downloads. 

A logical first step in 
developing new phenomenological 
information and predictive capabil- 
ity for this class of rotor-body 
interference aerodynamics would be 
to study the two-dimensional sec- 
tion characteristics cf a wing in 
the wake ot a rotor, or even more 
simply, to study an airfoil placed 
normal to an oncoming uniform flow 
(Fig. 3). This novel configuration 
is the basis of the present com- 
bined theoretical and experimental 
investigation. A special wind 
tunnel experiment has been per- 
formed for this problem, and two 
new modern computational aerody- 
namics methods have been explored 
to complement the measurements. 

The resultant two-dimensional data 
are expected to approximate the 
section characteristics that are 
being measured in a separate inves- 
tigation at various spanwise sta- 
tions of the actual XV-15 Tilt 
Rotor Aircraft in hovering flight. 


(B) TWO-DIMENSIONAL AIRFOIL 

Fig. 3. The rotor download problem 
and a strip-theory approximation. 


SECTION A-A 

An important aspect of the 
present investigation is the com- 
bination of experimental information 
and computational analysis. The 
experiment provides some definitive 
facts about the real separated vis- 
cous flow, but it has specific 
limitations with respect to wind- 
tunnel wall corrections, Reynolds 
number, and the limited number of quantities that are feasible to measure. 
These limitations are easily overcome by the numerical methods, and in 
addition, the effort involved in changing the computer input to modify the 
airfoil shape is much less than that of modifying physical wind-tunnel 
models. The flexibility to change the geometry at will and to examine the 
flow-field solutions in detail leads to a better theoretical understanding 
of the physics of the problem. However, the physical modeling and approxi- 
mations of the numerical methods have to be examined and verified, and 
further improvements are required to determine the absolute values of the 
air ^ads with confidence. We shall show that while both the experimental 
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and numerical approaches have definite shortcomings, the proper combination 
of computations and measurements gives more information than can be derived 
from either method alone. 


II. Experimental Investigation 

Two-dimensional models of the XV-15 wing with various flap and 
leading-edge configu. stions were tested between endplates in the U.S. Army 
Aeromechanics Laboratory 2x3 meter subsonic wind tunnel. Figure 4 shows 
the wing sections that were used to obtain the results given in this 
paper. The airfoil shown (a modified NACA 64A223 profile with a 25% plain 
flap) represents the XV-15 aircraft; 30% and 35% trailing-edge flaps were 
also tasted. The trailing-edge tlaps were deflected in 15° increments up 
to 90°. The modified leading edge was designed on the basis of preliminary 
calculations which revealed that the drag characteristics are highly sensi- 
tive to the surface curvature distributions in certain critical regions on 
the upper, or "windward," side of the airfoil. 




MODIFIED 
LEADING EDGE 



Fig. 4. Sketch of the models tested. 

The chord of the basic model., (with <Sp * 0 and no leading-edge 
modifications) was 0.31 m. This represented a difficult compromise 
between the requirements to maximize Reynolds numbers, minimize wind-tunnel 
blockage and wall effects, and minimize three-dimensional effects. The 
Reynolds number for the airfoil results presented herein was 10 6 . The 
results for a number of additional configurations and for ranges of angles 
of attack and Reynolds numbers are given in Ref. [2], along with further 
details of the experimental setup. 

In addition to the airfoils, two wedge-shaped models having 
equilateral-triangle cross sections, with c ■ 0.22 and 0.31 m, were tested 
with the apex pointed both forward and rearward. These two orientations 
produced values of Cq that were comparable to or greater than the values 
for the various flap settings. The data from these models were essential 
in developing and validating the test techniques, as explained in Sec- 
tion IV. The Reynolds number based on c varied from 0.4 to 1.3 x 10 6 for 
the triangles. 


ORIGINAL PAGE fS 
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A typical model installation is shown in Fig. 5. The spar of the 
models was cantilevered from the frame of a force-and-moment balance 
beneath the floor of the test section. A turntable in the balance frame 
allowed the model angle of attack to be adjusted ±20° from perpendicular to 
the free stream. Large endplates, based on the observations and recommen- 
dations of Ref. [3], were installed 0,31 m from the wind-tunnel floor and 
ceiling to minimize the interaction between the tunnel-wall boundary layer 
and the wake of the model. The small gaps between the ends of the models 
and the endplates were not sealed. 


TOP VIEW 



AIRFLOW 


SIDE VIEW 



Fig. 5. Wind-tunnel installation. 


The average aerodynamic forces were derived from the wind-tunnel 
balance system and static pressure distributions were measured at three 
spanwise locations on the airfoils. For the triangles, only two static 
pressure taps were installed on each face at each of the three spanwise 
locations. Other measurements included oil flow and wool tuft visualiza- 
tions of the separation patterns and wake surveys with fast-response pres- 
sure transducers. 


These flow visualization studies and spanwise traverses, along with 
the measured spanwise pressure distributions, indicated thac: the flow was 
uniform in the vertical direction to within the accuracy of the measure- 
ments even though the aspect ratios of the large triangle and the airfoil 
models were only 5. Fuithermore, no evidence was found of spanwise cells 
in the wake structure. The rather large blockage ratio (up to 10% based on 
the frontal area between the endplates) was a matter of concern and special 
study, but as indicated in Section IV, the corrected data for the triangles 
agreed well with previously published results. 
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III. Numerical Methods of Analysis 


ORIGIN/,;. . 

OF POOR QUALslY 

The first aerodynamic code to be considered was developed by the 
second author (Refs. 14,5]) to calculate the separated flow, wake, and 
fluctuating airloads on two-dimensional bluff bodies or airfoils at arbi- 
trary angles of attack and at high Reynolds number. This technique com- 
bines an integral boundary layer with a discrete vortex method for the 
outer flow. The two-dimensional vorticity-conservatlon equation is solved 
in a Lagrangian formulation, wherein the vorticity field is represented as 
the sum of local patches, or "blobs," of vorticity which retain their 
strength in time and are convected by the flow. The vortices are intro- 
duced along the walls of the body at each time step, and their positions 
at subsequent time steps are determined numerically using the Adams- 
Bashforth-2 multistep time-integration scheme. The resultant velocity 
field computed by the Biot-Savart law is used in a boundary-layer calcula- 
tion to determine the separation point. Viscous diffusion is neglected 
outside the attached boundary layer and no empiricism is introduced. 

The present discrete-vortex method does not use conformal mapping; 
hence, it can treat arbitrary shapes and multiple bodies. This feature has 
been exploited to include flat surface.® that represent wind-tunnel walls, 
which is essential for quantitative comparison with the experiment. The 
code gives the complete time-dependent development of the entire flow field, 
including vortex shedding. However, it requires relatively large computa- 
tional resources for large numbers 
of vortices; a typical case 
requires 10-15 min CPU time on the 
Ames Cray IS computer. Figure 6 
shows a typical result. Here the 
dots are the individual vortices 
and the contour lines are the 
instantaneous streamlines. 

The second code is a recent 
adaptation by the fifth author of 
the program VSAERO (Ref. [6]), 
which combines an efficient three- 
dimensional, unsteady potential- 
flow panel method with a free- 
streamline representation of the 
separated zone. Planar quadri- 
lateral panels are used to repre- 
sent the body and wake surfaces. 

Each panel has a constant source 
and doublet distribution and a 
central control point where an internal Dirichlet boundary condition is 
applied. Large regions of separated flow are modeled in the manner of 
the CLMAX program (Ref. [7]), which assumes an inviscid wake with total 
pressure that is less than the free stream value. The separated-wake 
region is enclosed by a pair of constant-strength vortex sheets. 

The calculations for this method proceed as follows: an initial 

solution is assumed, including the separation points and the shape of the 
wake, and then the solution is stepped forward in time. The dividing 
streamlines between the potential-flow zones and the wake region are trans- 
ported with the local outer flow at each time step using the calculated 
velocities of points on the wake surfaces, thereby satisfying the condition 
that the wake vortex sheets be force-free. In the present case, no attempt 



Fig. 6. The instantaneous flow field 
computed around an NACA 4421 airfoil. 
The arrow indicates the position, 
magnitude, and direction of the 
resultanc force vector. 
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is made to model the unsteady vortex shedding phenomenon; rather, the code 
is run until the solution converges to a steady-state solution that approx* 
imates the time-averaged separated flow around the body, usually within 
10 time steps. 

By its nature, this code contains more approximations and empiri- 


cism than the discrete-vortex approach 



Fig. 7. Pressure distribution on the 
NACA 64A223-M airfoil at a = -90° 
and 6p = 60°. 


However, the two-dimensional ver- 
sion requires only about 100 sec of 
CPU time per case on a Prime 550 
computer to converge to a steady 
solution; this would be equivalent 
to one or two seconds on the 
Cray IS. Figure 7 shows a typical 
result. The large suction peak on 
the upper surface at x/c * 0.85 
is due to the rapid expansion 
around the shoulder of the flap. 

Both of the codes used in 
this study have been adapted to 
include airfoils with flaps, 
multiple-element airfoils, and 
exterior boundaries such as wind- 
tunnel walls. This latter capabil- 
ity is essential for detailed com- 
parison with the experiment and for 
verifying the wall corrections that 
were applied to the data. 


IV. Results and Discussion ~ Triangles 

A. Validation of the Experiment, Including Blockage Corrections 


When a model is tested in a closed-section wind tunnel, it creates 
a blockage that accelerates the local flow and increases the drag. These 
effects are known to be proportional to the drag and the physical size of 
the model. For two-dimensional tests, Allen and Vincent (Ref. [8]), 
Pankhurst and Holder (Ref. [9]), and Maskell (Ref. [10]) give theoretical 
blockage corrections that take the following form for bluff bodies. 



( 1 ) 


C 

c 

Po 


1 

T 



( 2 ) 


where e is a constant, b is the ratio of the lateral dimension of the 
model to the lateral dimension of the wind tunnel, C D is the two- 
dimensional drag coefficient in free air, and is the measured uncor- 

rected value of the drag coefficient. Also, C p and C p& are the corrected 
and uncorrected values of the pressure coefficient, respectively. 
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References [8-10] give theoretical values of e ranging from 0.50 
to 0.96, and numerous previous wind-tunnel measurements on cylinders of 
various cross sections suggest values between these limits. Further sup- 
port for this approach was obtained from the present discrete-vortex compu- 
tational method. The results for airfoils at a = -90° between solid 
walls were found to correlate well with Eq. 1 for blockage ratios in the 
range 0 < b < 0,20, giving c = 0.65 ±0.05. 

Therefore, the form of Eq. 1 seems to be appropriate here, provided 
a reasonable estimate of e can be obtained. This empirical constant was 
obtained from the present data for the two different sizes of triangles 
tested at the two different orientations, giving four combinations of the 
product bC D . The corresponding free-air drag coefficients have been well 

documented in Ref. [11]; C D * 2 00 for the blunt face forward and 
Cj) = 1.30 for the apex forward. The measured values for 15 combinations 
of Reynolds numbers, sizes, and orientations then yielded e * 0.596, with 
a standard deviation of only ±0.024. As this value is in good agreement 
with the various independent studies cited above, it was used to correct 
the airfoil data described in the following section. 

The experimental results for the triangles are summarized in Table 1. 
Only the average values of the various experimental quantities are listed, 
as they were found to be independent of Reynolds number, to within the 
experimental uncertainty. All of the results a**e in excellent agreement 
with Hoemer (Ref. [11]), with the exception that the Strouhal frequency 
for the triangles with the apex forward (St - fb/U^) does not correlate 
with his empirical formula using Cj). It is interesting to note that the 
corrected base pressure coefficient, Cp^* is essentially independent of the 

orientation of the triangle, even though Cp and St are not. 


TABLE 1. SUMMARY OF RESULTS FOR THE EQUILATERAL TRIANGLES 


Configuration 

C °0 

C D 

C D 

(Ref. 11) 

^base 

Cp 

base 

(Ref. 11) 

St 

St 

(Ref. 11) 


2.31 

1.99 

2.00 
1 . 98 1 

-1.17 

-1.13 1 

0.128 

** 

CM 

d 

- ► 

2.25 

2.03 

2.00 
1 . 98 l 

-1.24 

-1.13 1 

0.123 

0.123 2 


1.41 

1.29 

1.30 

-1.12 

-1.13 1 

0.200 

0.173 2 

- ◄ 

1.38 

1.30 

1.30 

-1.18 

-1.13 1 

0.204 

0.173 2 


x Flat plate normal to flow. 
2 Hoemer: St - 0.21 Cg 3 '**. 
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ORIGINAL F.V.£ r H 

B. Validation of the Computational Methods OF POOF •• * 


The data in Table 1 can now be used to evaluate the present compu- 
tational methods. Figure 8 shows the free-air drag coefficient for tri- 


angles as a function of the semi- 
vertex angle. As seen from the 
figure, the simpler panel method 
gives reasonably good results, 
but the vortex method does not. 

On the other hand, the computed 
Strouhal shedding frequencies of 
the vortex shedding (not shown) 
were within a few percent of the 
experimental values. The main 
difficulty with the vortex method 
seems to be that the computed 
base pressure is much too large, 
Cp^ = -2.0 vs -1.2 in the experi- 
ment, and this is responsible for 
the excessive values of Cp. The 
panel method gives approximately 
the correct base pressure for 
both of the triangles tested. 


HOER^ER 

• EXPE .WENT (CORRECTED) 
A VORTEX CALC 
□ PANEL CALC 



■ ■ i . —i 1 

0 30 60 90 

SEMI-VERTEX ANGLE. >1 


Fig. 8. Drag coefficients for tri- 
angles as a function of semivertex 
angle. 


Despite efforts to deter- 
mine the deficiency in the vortex 
method, the reasons for it remain 
unknown. The method was shown in 
Refs. [4] and [5] to predict 
dynamic stall on an oscillating 

airfoil and the flow field of a circular cylinder reasonably well, but was 
less successful in determining the drag of a square cylinder, which was 
found to be too low. A sensitivity study of the numerical parameters such 
as time step, number of points used to define the body, number of vortices, 
and vortex core radius has thus far failed to reveal any clear trends. 


However, it should be mentioned that previous Investigators of vor- 
tex methods have found It necessary to reduce empirically the circulation 
of the vortices after they leave the body (cf. Sarpkaya [12]). This la 
often argued as modeling vortlcity dissipation due to viscosity even though 
vortlcity inly diffuses within the framework of the Navier-S tokes equations. 
To test the importance of diffusion, the effects of viscosity were simu- 
lated in test calculations by means of Chorin's "random walk" (Ref. [13]). 
Changes in the base pressure were found; but only by simulating low values 
of Reynolds number on the order of 100, was the drag reduced to approxi- 
mately the experimental values. The remaining possibility is that the 
vortlcity in the wake becomes highly three-dimensional and this may somehow 
reduce its effective induced-velocity field in the plane of the mean flow. 

At present, however, the present vortex method cannot be considered relia- 
ble for quantitative predictions without empiricism, although it may be 
valuable for predicting trends. 


V. Results and Discussion — Airfoils 


The primary objectives of this investigation wore to determine the 
aerodynamic characteristics of the XV- 15 airfoil section with various flap 
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O BASED ON FULL CHORD 

□ BASED ON FRONTAL AREA 


kf, deg 

Fig. 9. Measured drag coefficients 
of the airfoil with 25% trailing- 
edge flap as a function of flap 
deflection angle. 


deflections and to determine whether 
the tilt-rotor downloads could be 
reduced by improved airfoil and flap 
designs. Figure 9 shows the mea- 
sured results for the drag coeffi- 
cient for the XV-15 profile. These 
resales are based on two different 
reference areas; (P ea of 

the basic airfoil v5‘h no 1 ' 

deflection, and (2) he actu^ ' pro- 
jected frontal area. In the latter 
case, the purely geometrical effect 
of reducing the surface .»*ea normal 
to the flow has been eliminated, and 
the variation in drag coefficient 
defined in this way is due to the 
modified aerodynamics alone. The 
shapes of the faired curves between 
6p = 45° and 60° were determined 
from crossplots of the results at 
other angles of attack. 


The results show that for flap deflection angles less than 60°, the 
total drag decreases significantly more than could be explained on the 
basis of the reduction in frontal area. However, for 6p > 60°, the total 
drag remains approximately constant, and based on frontal area 

actually increases. Tuft and oil-flow visualization revealed this to be a 
result of flow separation occurring on the flap just downstream of the 
shoulder of the flap and ahead of the trailing edge. This produced a wider 
wake behind the airfoil, less suction on the front face of the airfoil, a 
somewhat lower base pressure, and higher drag. 



. 10. Measured drag coefficients 
several airfoil configurations: 
is based on the chord of the 
airfoil with 6p - 0, xh is the 
location of the flap hinge axis* 


Fig 

for 

CD 


Figure 10 shows a comparison 
of the experimental results for all 
three flap sizes and for the modi' 
fied leading edge, as depicted in 
Fig. 4. The larger flaps supported 
attached flows to larger flap- 
deflection angles, with correspond- 
ingly lower values of Cp. A mini- 
mum value of Cp • 1,0 is estimated 
for the unmodified leading edge. 
However, the modified leading-edge 
reduced the minimum drag coefficient 
to only 0.64 with 6p ■ 60 s . These 
results indicate ••he potential value 
of wing modifications in reducing 
the downloads on the tilt-rotor 
aircraft . 


Figures 11 and 12 show the 
computations for the airfoil with a 
25% traillng-edge flap in comparison 
with the measured data from Fig. 9. 

It is clear from Fig. 11 that the 

quantitative prelections leave something to be desired, especially regarding 
the results of the vortex method. However, Fig. 12 indicates that the 
trends can be predicted quite well; for the drag coefficient ratio, the 
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EXPERIMENT 

NACA 64A223M 

O 5 
a 

A 

VORTEX CALC 

NACA 4421 

V 

VORTEX CALC 

NACA 64A223M 


□ 

PANEL CALC 

NACA 64A223M 



FLAGGED SYMBOLS - SEPARATION 
FIXED AT SHOULDER OF FLAP 


• EXPERIMENT 
V A VORTEX CALC 
□ PANEL CALC 


Fig. 11. Comparison of measured 
and calculated drag coefficient 
as function of flap deflection 
angle. 


Fig. 12. Measured and calculated 
drag coefficients normalized by the 
values for 6p ■ 0. 


panel method agrees with the mea- 
surements to within experimental 
accuracy. It is also interesting 
to note the following ratios of 
Cd for the airfoil with the 
modified leading edge and 
6 F » 60° compared to the basic 


airfoil with 


experi- 


EXPERIMENT, UPPER 
EXPERIMENT, LOWER 
SEPARATION ON FLAP 
SEPARATION AT T. E. 


ment, 0.37; vortex method, 0.53; 
panel method, 0.49. 


Figure 13 Lhows a compari- » JJ 

son of the measured pressure dls- ^ _ II 

tribution vs that predicted by X 

the panel method. In its present II 

preliminary form, the unsteady t l , * 

panel code does not include a 0 t/e 6 10 

boundary-layer calculation and 

the separation point must be pre- __ . ... 

scribed. However, it is clear Fig * l3 ‘ Measured and predicted 

from the pressure distributions TT™* di,trlbutions for 

that the flow would not remain f 

attached all the way to the 

trailing edge in this case; therefore, fixing the separation point on the 
shoulder of the flap gives somewhat better results. As the location of the 
separation point is like’y to be less obvious in other esses, the logical next 
step in the development of the method is to include s boundary- layer model. 
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VI. Co nc 1 ud i ng Remarks 


This investigation has produced new insight and quantitative infor- 
mation about the airloads on bluff bodies, with particular relevance to 
the complicated aerodynamic interference betwe n the rotors and the wing 
of tilt -rotor aircraft configurations. Both the calculations and the 
experiment show that the drag of an airfoil normal to the oncoming flow 
decreases as the flap deflection angle increases, up to the point where the 
flow begins to separate on the flap. Furthermore, the reduction in drag is 
considerably more than would be due merely to the reduction in the pro- 
jected area normal to the flow. However, the drag increases with increas- 
ing 5 y once separation appears on the flap, and this occurs well before 
the upper surface of the flap is aligned parallel with the free-stieam 
f low. 


The results for the airfoil model of the XV- 1 5 wing (Fig. 9) help to 
explain the behavior of the flight-test data in Fig. 2, as discussed in the 
Introduction. The hover performance is, of course, directly affected by 
the download on the wing (that is, by C D ), which depends strongly upon 
fip- Therefore, it is clear that the gross-weight capability in hover 
should increase with increasing flap deflection until flow separation 
begins tu occur on the flaps, but that excessive flap deflections would 
increase the wing download and decrease the maximum gross weight in hover. 
The dashed line in Fig. 2 is the original fairing of the data as presented 
in Ref. [1), whereas the solid line reprei* 'ts the revised estimate based 
on the results of the present investigation. 

The combined results of this investigation also indicate that sig- 
nificant further improvements in hover efficiency could be attained by 
careful design of the wing sections. The most obvious possibilities 
include the use of larger flaps, flaps with larger radii of curvature at 
the shoulder, and appropriate changes in the curvature distribution in the 
leading-edge region of the wing. Multielement airfoils were not considered 
in this paper, but some further drag reductions due to extra devices are 
described in Ref. [2]. 

The experimental techniques that were developed with the aid of the 
computational methods and the triangle models appear to have been very suc- 
cessful in dealing with and correcting for the relatively large blockage 
ratios of the airfoils. As a side benefit, some additional data have been 
added to what exist in the general literature for arbitrary bluff bodies. 
The quantitative accuracy of the discrete-vortex computational method 
turned out rather disappointing, although it was still useful. The reasons 
for its failure and a means of introducing suitable empiricism should be 
examined further. On the other hand, the panel method with the free- 
streamlire representation of the separated wake is quite promising. As 
mentioned in Section V, some means of estimating the separation points by 
boundary-layer theory ia essential if the method is to be used to design 
the optimum curvature distributions in the leading-edge region and near the 
shouLder of the flap. Furthet efinemente are needed to enhance its quan- 
titative accuracy as well, but the method seems to be a good. Inexpensive 
engineering tool to study complex flow problems. 

Finally, it should be emphasised that the piesent investigation was 
concerned entirely with a two-dimensional approximation to a very compli- 
cated tnree-dimensional aerodynamic-interference problem. The results 
Indicate that considerable improvements are possible anF practical in the 
area of rotor-induced downloads, but similar studies for realistic 
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rotor-wing combinations will be essential to help the tilt-rctor concept 
achieve its full potent { ?1. 
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